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ABSTRACT 



We analyze a sample of galaxies chosen to have F2A^m. > 0.5mJy and satisfy 
a certain IRAC color criterion. IRS spectra yield redshifts, spectral types, and 
PAH luminosities, to which we add broadband photometry from optical through 
IRAC wavelengths, MIPS from 24-160 /xm, 1.1 millimeter, and radio at 1.4 GHz. 
Stellar population modeling and IRS spectra together demonstrate that the dou- 
ble criteria used to select this sample have efficiently isolated massive star-forming 
galaxies at z ^ 1.9. This is the first starburst-dominated ULIRG sample at high 
redshift with total infrared luminosity measured directly from FIR and millimeter 
photometry, and as such gives us the first accurate view of broadband SEDs for 
starburst galaxies at extremely high luminosity and at all wavelengths. Similar 
broadband data are assembled for three other galaxy samples - local starburst 
galaxies, local AGN/ULIRGS, and a second 24/im-luminous 2r ~ 2 sample dom- 
inated by AGN. Lpah/Lir for the new 2; ~ 2 starburst sample is the highest 
ever seen, some three times higher than in local starbursts, whereas in AGNs 
this ratio is depressed below the starburst trend, often severely. Several pieces 
of evidence imply that AGNs exist in this starburst dominated sample, except 
two of which even host very strong AGN, while they still have very strong PAH 
emission. The ACS images show most objects have very extended morphologies 
in the rest-frame UV band, thus extended distribution of PAH molecules. Such 
an extended distribution prevents further destruction PAH molecules by central 
AGNs. We conclude that objects in this sample are ULIRGs powered mainly by 
starburst; and the total infrared luminosity density contributed by this type of 
objects is 0.9 - 2.6 x WL^/Mpc^. 

Subject headings: cosmology: observations — galaxies: dust emission — galaxies: 
mid-infrared 



INTRODUCTION 



Formation of massive galaxies provides a critical test of theories of galaxy formation 
and evolution. Before modern deep observations, the most massive galaxies known were 
local elliptical galaxies \y ith no ongoing star formation. The classical model for these objects 
(e.g., lEggen et al.l Il962l ) was monolithic formation at high redshifts, followed by passive 



evolution. A more recent galaxy formation theory in the Cold Da rk Matter (CDM) paradigm 
( IBlumenthal et al.lll984 IWhite fc Frenklll99ll : ICole et al.ll2000l ) predicts quite the opposite 
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scenario: a galaxy-galaxy merging-tree model. In this scenario, small galaxies formed early in 
cosmic time, and massive galaxies were assembled later at much lower redshifts by a series of 
mergers. Observations of local Ultra-Luminous InfraRed Galaxies (ULIRGs, > 10^^ Lq 



detected by IRAS are consistent with the merger theory. Most loca . 



morphologies, consistent with being mergin g systems (IKim et al. 



stages of me r ging have r ' light profiles (1 James et al 



1999 



■ ULI RGs have disturbed 
2002h. ULIRGs in later 



Rothberg fc JosephI \20Q4 \. 



Genzel et al.l ( l200ll ) and iTacconi et al.l (120021 ) measured local ULIRG dynamical masses and 
found an average of 10^^ Mq. These features are consistent with numerical simulation studies 
of galaxy mergers, indicating that lo cal ULIRGs are merging sys t ems transforming gas-rich 



galaxies into e l liptical galaxies (jKormendy &: Sanders! Il992l : iMihos fc HernquistI Il996 
Barnes et aDll996l : IXacconi et aDl2002[ ). 



The story is different at z > 2. Deep near-infrared surveys ( Franx et al.ll2003l : iGlazebrook et al 
20041 : iMcCarthv eTaPbooi IPaddi et allboosl : iLabbe et al.lboosh have identified apparently 
luminous passive galaxies already in place at 2; ~ 2, implying that they formed at even higher 
redshifts. The existence of galaxies with M^, > lO'^^ M© at high reds hifts may challenge the 
merger theory of forming such objects at lower redshifts. However, ICattaneo et al.l (l2008l ) 
used a semi-analytic model to show that significant numbers of > IO^^Mq galaxies were 
in places by 2; ~ 2 but many also formed at lower redshifts, that is, there was a whole 
" downsizing" trend for massive galaxies to form their stars early, but it is merely statistical 
not absolute. Possibly consistent with this is the fact that the contr ibution of LIRGs and 
ULIRGs to the total IR luminosity density is more than 70% at z_= 1 (ILe Floc'h et al.ll2005[) 



compared to a negligible percentage locally (ISanders fc Mirabellll996l: iHuang et al 



2007 



More over, the redshift surveys for Sub-Millimeter Galaxies (SMGs) by lChapman et al.l (12003 



2OO5I ) reveal a rapidly evolving ULIRG population at 1.7 < z < 2.8. Such strong evolution is 
also seen in ULIRGs selected with BzK color and MIPS 24 /xm fiux at 2; ~ 2 with their num- 
ber density apparently 3 orders of magnitude higher than the local number density. Thus 
local ULIRGs would well be the tail end of earlier intense activity. 

The Spitzer MIPS 24 /xm ba nd has been yery effective in probing infrared emission frona 



galaxies at redshifts up t o z ~ 3 JHuang et ahlboosl : I Webb et al.l 



Rigopoulou et al. 



2006 



Daddi et al. 



2005 



2007a 



2006 



d). iPapovich et al 



Papovich et al. 



2006 



(120061 ). IWebb et al. 



(l2006h. lDaddi e t al. N2007J bl)7 and lDey et al. I J2008h argued that 24 /im emission from galax- 



ies at 2 < z < 3 is powered by both active galactic nuclei (AGN) and star formation. Spec- 
troscopic o bservations of a f e w 24 nm lurainous SMGs and Lyman break galaxies (^LBGs) at 



1 < z <3 (ILutz et al.ll2005l : iHuang et al.ll2007al : iLe Floc'h et all 120071 : IVahante et al.ll2006 



^IR 



Js /zm 



L^df (jSanders k Mirabe]||l996l ) 
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Pope et al.l 120081 ) with the Infrared Spectrograph (IRS) on Spitzer support this view, show- 
ing both strong continua and emission features of Polycyclic Aromatic Hydrocarbons (PAH) 
in the rest-frame 6 < A < 10 fim. Systematic infrared spectroscopic surveys of 24 /xm lu- 
minous but optically faint sources (IHouck et al.ll2005l : lYan et al.ll2005l ) reveal a dusty, 2; ~ 2 
AGN population not detected in optical surveys. Most of these AGNs are ULIRGs with 
power-law spect ral energy distributi on s (SEDs) in the mid-infr ared and deep silicate absorp- 
tion at 9.7 /im ( ISajina et al.ll2007al ). IWeedman et al.l (l2006bl ) observed a sample of X-ray 



AGN with similar properties though genera lly less silicate absorpti on. Optically-faint radio 
sources are a mix of AGN and starbursts (IWeedman et al.l l2006al ) but are predominantly 



AGN. In general, optically-faint obje cts have weak infrared spectral emission features, and 
most objects are likely to be AGN (IWeedman et al.l l2006d ). Ho wever, not all 24 um lu- 
min ous objects a t z ~ 2 — 3 are AGN-dominated. For example, IWeedman et al.l (l2006bl ) 



and iFarrah et al.l (120081 ) also identified samples with an apparent 1.6 /im stellar peak in the 
IRAC 4.5 or 5.8 /xm bands. Both samples show a very narrow redshift distribution due to 
a selection of the MIPS 24 /zm band to^ yard strong 7. 7jU.ni PAH emission at z ^ 1.9. IRS 
spectroscopy of 24 /im-luminous SMGs (ILutz et al.l l2005l: iMenendez-Delmestre et al.l 120071 : 
Valiante et al.ll2006l : iPope et al.ll2008l : iMenendez-Delmestre et al.ll2008l ) shows similar spec- 
tral features, namely strong PAH emissio n in objects with a 1.6/im stellar emission bump 
( IWeedman et al.ll2006bl : iFarrah et al.ll2008l ). indicating intensive star formation in both types 
of objects. 

This paper presents an IRS spectroscopic and multi-wavelength study of a ULIRG sam- 
ple at z ~ 2. The sample conies fro m the All- wavelength Extended Groth-strip International 
Survey (AEGIS, iDavis et al.l 120071 ). which consists of deep surveys ranging from X-ray to 
radio wavelengths. Selection of our sample catches a starburst-dominated phase of ULIRG 
with Lj^ > 5 X IO^^Lq, which is very rare among local ULIRGs. In this paper, we will study 
their properties including star formation, stellar masses, AGN fractions, and contribution 
to the universe's star formation history. §2 describes the sample selection. The IRS spec- 
troscopic results are presented in §3, and §4 contains an analysis of stellar populations, star 
formation rate, and AGN fraction. §5 summarizes our results. All magnitudes are in the AB 
magnitude system unless stated otherwise, and notation such as "[3.6]" means the AB magni- 
tude at wavelength 3.6 fim. The adopted cosmology parameters are Hq = 70 km s"^ Mpc~^, 
Qm = 0.3, = 0.7. 
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SAMPLE SELECTION 



We wish to study the mult i- wavelength properties of star-forming galaxies at z ~ 2. 
There are many ways of using optical and NIR colors to select such a sample. The sam- 
ples with the best spectroscopic confirmation are the UgR color-selected BM/BX sources 



( Steidel et al.l 120041) . which have estima t ed ty pical stellar masses of about 10^ ~ 10^° M( 



( IShapley et al.l l2005l : iReddv et al 



2006 







20081 ). The average 24 fim flux density for these 



sources is 42 ± 6 fiJj (IReddy et al.ll2006l ). which suggests modest rest frame mid-IR lu 
minosities, consistent with LIRGs (-L/_r < 10^^ -^o)- A different sample, ba sed on near 
infra r ed color selec t ion, i s the Distant Red Galaxies (DRG, (J — K)yega > 2.3) (IFranx et al. 
20031 : iLabbe et al.l l2005l ). These galaxies are redder and dustier tha n the UV-selected 



BM/BX sources and are believed to be more massive than W^^Mq flLabbe et al.l l2005l 
Papovich et al.l 120061 ). Dusty DRGs have estimated total iii frared luminosity in the range 
10^1 <LiR< W^Lq dWebb et ahlbood : IPapovich et al.lbooel ). A third sample JDaddi et al. 



20051) uses B z K col ois to select galaxies at 2; ~ 2; massive galaxies are mid-IR luminous. 
Reddy et al.l (120051 ) compared BM/BX, DRGs, and BzK gal&ides and found that BzK galax- 
ies include rnost D RGs and BM/BX galaxies. This comparison is nicely shown in Fig. 9 of 
Reddv et al.l J2005h . 



An independent way to select galaxies at 2; > 1.5 is to use IRAC colors. In 1.5 < z < 3.0 
the four IRAC bands probe the rest-frame NIR bands where galaxy SEDs have similar shape 



thus the IRAC colors are very robust in determining redshift in this range (iHuang et al.ll2004 



Papovich et al.l 120081 ). The 1.6 /xm stellar emission bump can be used to separate galaxies 
ai z = 1.5. At 2; < 1.5, the IRAC 3.6 and 4.5 /xm bands sample galaxy SED from the Jeans 
tail of cold stars, thus the [3.6] — [4.5] < 0. At z > 1.5, the 1.6 fim stellar emission bump 
begins to move into the IRAC 4.5 n m band, making the [3.6] — [4.5] > 0. The color criteria 
is set based on the M82 SED model JHuang et al.ll2004h 



0.05 < [3.6] - [4.5] < 0.4, and 



- 0.7 < [3.6] - [8.0] < 0.5 (2) 

Both color ranges are corresponding to redshift range of 1.5 ^ z < 3.3. The red color 
cut in both equations rejects power-law AGNs and star forming galaxies at 2; > 3.3 whose 
7.7/im PAH shifts out of the IRS wavelength range. The selection is based on the rest- 
frame NIR colors, and is thus less affected by dust extinction, stellar ages, and metallicities. 
Figure [T] compares IRAC two-color plots for the DEEP2 galaxies with z < 1.5 to LBGs, 
DRGs, and BzK galaxies. The color criteria of equations 1 and 2 include most galaxies in 
the 1.5 < z < 3.0 range. Although 8 /xm detection is required for this selection method, 
selecting at this wavelength has some additional advantages. Chosen at roughly 2-3 fim 
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rest-frame, the sample selection is immune to dust reddening and is roughly equivalent to 
a stellar- mass-selected sample. It is thus ideal for stud ying luminous, potentially massive 



galaxies flHuang et all 12004 l2005l : IConsehce et al.l 120071 ) 



The specific IRS targets for this program were selected from a 24 /xm sample (jPapovich et al 



20041 ) in the EGS region. These sources show a clump at the predicted colors for 1.5 < 2; < 3 
in Figure [21 but redshifts were not known in advance. Individual targets were selected to 
have IRAC colors satisfying Eqs. 1 and 2 and also -^24^111 > 0.5 m Jy. In addition, e ach 
candidate was visually examined in the deep Subaru R-band image (lAshby et al.ll2008l ) to 
avoid confused or blended targets. With these criteria, 12 targets were selected in the 2°xl0' 
EGS region for IRS observation. Table 1 lists the sample galaxies. In this redshift range, 
most sources will be either ULIRGs with total infrared luminosity Ljji > 10^^ Lq |^ or AGNs 
with high mid-IR luminosities. For convenience, the galaxy nicknames used in the Spitzer 
database are used in this paper, but these do not follow proper naming conventions and 
should not be used as sole object identifiers. Proper names are also given in Table 1. 

Most of the previous IRS surveys of IR luminous sources at ^ ~ 2 have used rather 
different selection criteria. Table 2 summa r izes t he sample criteria for various other IRS 
surveys. iHouck et al.l (120051 ) and lYan et al.l (120051 ) used extreme optical-to-24 fim color to 
select dusty objects. Objects in these samples have much redder [3.6]- [8.0] IRAC colors 
than the majority of 24 fim sources (Fig. [2]) and are mostly A GNs as shown by their strong 



powe r-l aw continua, but weak o r absent PAH emission features ( IHouck et al.ll2005l : lYan et al. 



2OO5I ). IWeedman et al.l (l2006bl ) selected AGN using similar criteria. They also selected a 
separate starburst-dominated sample at z ^ 2 based on the stellar 1.6 /im emission bump. 
The exact criterion required the peak flux density to be at either 4.5 /im or 5.8 /im, thus 
rejects low-redshift galaxies and AGN with strong power-law SEDs . The re sulting sample is 



very similar to ours though overall a bit redder (Fig. [2]). All objects in the IWeedman et al. 
starburst sample show strong PAH emission features. 



^ The total infrared luminosites (L ib) for our sample are calculated with MIPS 24, 70, 160 /xm and 
l.lmm flux densities and Chary-Elbaz |Charv fc ElbajboOlllPaddi et al.ll2007al lbl) SED models. Details are 



given m 
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3. IRS Spectroscopy 



3.1. Observations and Data Reduction 



IRS observations of this sample are part of the GTO program for the Spitzer/YRAC 
instrument team (program ID: 30327). Objects were observed only with the IRS Long- 
Slit Low- Resolution first order (LLl) mode, giving wavelength coverage 20 < A < 38 /im 
with spectral resolution 60 < A/AA < 120. The wavelength coverage corresponds to 6 < 
A < 13 /im in the rest-frame for galaxies aX z ^ 2. This wavelength range includes strong 
PAH emission features at 7.7, 8.6, and 11.3 /im and silicate absorption from 8 to 13 /im 
(peaking near 9.7 /im). Detecting these features permits redshift measurement and study of 
dust properties. Total exposure tirn e for each target was based on its 24 /xm flux density. 
Mapping mode (ITeplitz et al.l 120071 ) was used to place each object at 6 positions spaced 
24" apart along the 168" IRS slit. This not only gives more uniform spectra for the target 
objects, rejecting cosmic rays and bad pixels, but also increases sky coverage for possible 
serendipitous objects around each target. Table 1 gives the target list and other parameters 
for the observations. All data were processed with the Spitzer Science Center pipeline, 
version 13.0. Extrac tion of source spectra was done with both the SMART analysis package 
( iHigdon et al.ll2004l ) and our customized software. Lack of IRS coverage at A < 20 /^m for 
this sample is compensated with deep AKARI 15 /xm imaging (jim et al.ll2008l ). All objects 
except two outside the AKARI area are detected at 15 /zm providing measurement of the 
continua at rest-frame ~ 6 /im for galaxies ai z ^2. 

Figure [3]presents the IRS spectra. PAH emission features at 7.7 and 11.3 /xm and silicate 
absorption peaking at 9.7 ixm are detected from 10 sources, indicating a narrow redshift 
range of 1.6 < z < 2.1. The PAH emission features at 7.7 and 11.3 /xm show pronounced 
variations in their profiles and peak waveleng ths. Both 7.7 and 11.3 yum PAH features have 
at least two components ( jPeeters et al.l 120021 ). For example, the 7.7 ^m. PAH feature has a 
blue component at 7.6 /im and a red component at wavelength longwards of 7.7 /im. Thus 
different types of PAH spectral templates potentially yield different redshift measurements. 
To check this, we use two local MIR spectral templates with different PAH profiles, an average 
local starburst spectrum and an average local ULIRG spectrum to determine redshifts. Both 
templates yield very similar redshifts (Table 3). The starburst template fits all objects better 
with a typical 2% redshift uncertainty. EGS_b2 is identified at 2; = 1.59 with PAH emission 
features at 8.6 and 11.3 /im and the [Ne II] emission line at 12.81 /im. Redshift z = 2.03 
for EGS12 is confirmed by detecting Ha at 1.992 /im ( Figure HI) in a NIR spectrum taken 
with the MOIRC spectrograph on the Subaru telescope (lEgami et al.ll2008l ). The spectrum 
of EGS_b6, however, shows two emission lines at 27.7 and 31.1 fim that we are not able 
to identify consistently with any redshift. EGS_b6 is resolved to two objects 0"7 apart in 
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the HST ACS image (jPavis et al.l 120071 ). and an optical spectrum of this system shows two 
galaxies at z = 1.02 and z = 2.001. We therefore omit EGS_b6 from the sample for further 
analysis. The 24 fim images show several serendipitous objects in slits of all 12 targets, 
most of which are too faint to permit redshift identification. Only one source, EGS24a, have 
-^24 nm ~ 1 mJy. This object, found in the slit of EGS24, shows the silicate absorption 
feature at z = 2.12 (Fig. [3]). 

The redshift distribution of the sample (F i g. El) is very similar to that of the starburst- 
dominated ULIRGs studied by lWeedman et al.l (l2006bl ). e ven though ou r limiti ng flux density 
at 24/im is a factor of two fainter than theirs. Recently iFarrah et al.l (120071 ) use the same 
criteria to select a larger sample in the Lockman Hole region for the IRS observation and 
yield a very similar redshift distribution. The narrow distribution for starburst-dominated 
ULIRGs is due to the selection of strong 7.7 fim PA H emission by the MIP S 24 /xm band 
at z ~ 1.9. The peak of the redshift distributions for lWeedman et al.l (l2006bl ). iFarrah et al. 



(120071 ). and our sample is at this redshift, confirming the selection effect. On the other 
hand, lumi nous 24 /im sources with power-law SEP have a much wid er redshift range up 
to 2 ~ 3 JHouck et al.l l200,4 IVan et al.l boosi : IWeedman erall l2006bh . but t hev will not 
pass our IRAC colo r criteria or the "bump" SED criterion in IWeedman et al.l (l2006bl ) and 
Farrah et all J2007h . 



3.2. PAH Emission Features in ULIRGs 



The PAH features visible in the individual spectra of the sample galaxies are even 
more prominent in th e average spectrurn for the sample, as showed in Fig El which also 
stacks local starburst (IBrandl et al.l l2006l ) and ULIRG samples for comparison. The local 
ULIRG sample is divided int o Seyfert, LINER, and HII sub-samples according to their 
optical spectral classific ation (IVeilleux et al.lll999l ). PAH emission features are found have 
different feature profiles. iPeeters et al.l (|2002| ) classified profiles of each PAH emission feature, 
according to the peak wavelength, into 3 main classes: Class A, B, and C. PAH emission 
features are known to have more than one component in each feature. For example, the 
7.7 /xm PAH emission feature have two major components at 7.6 and 7.8 /im: Class A is 
defined as 7.6/im dominated PAH; Class B is the 7.8/xm component dominated PAH; and 
Class C is red component dominated with peak shifting beyond 7.8/xm. The 7.7 /im PAH 
in the local starburst spectrum appears to be more consistent with class A with the peak 
at wavelength shorter than 7.7 /xm. All local ULIRG spectra have a typical class B PAH 
profile, with a red wing extending beyond 8 /im. In §3.1, we already found that the starburst 
template fits each IRS spectrum of our sample better than the ULIRG template. It is not 
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surprising then that the average 7.7 /im PAH profile of our sample is more similar to the 
average starburst spectrum, thus consistent with the class A. Another significant difference 
is that our ULIRG sample has an average Ln 3y/,m /-^7 7/A]yil^ ra tio about twice as high as 



local ULI RGs but similar t o local starbursts (IBrandl et al.ll2006l ). We also plot the average 



spectra of lYan et al.l (120051) in Figure [6] for comparison. The average spectrum for strong 
PAH objects in lYan et al.l (120051 ) is more similar to the local Seyfert type ULIRG, implying a 
dominant AGN contribution in the spectra of their sample. We conclude from IRS stacking 
that the 7.7 /im PAH profiles and i^ii.3^m/-^7.7/xm ratios for the present sample are more 
consistent with those of local starburst galaxies rather than local ULIRGs. 

PAH emission features are a tracer of star formation, one of the energ y sources powering 



ULIRGs ( IGenzel et al.l Il998l : iRigopoulou et al.l Il999l : ISaiina et al.ll2007al ) . In o rder to sub- 



tract the local continuum, we adopted the method used by lSajina et al.l (j2007al ). fitting the 
5 < A < 15 /im spectrum into three components: the PAH emission features, a power-law 
continuum, and the silicate absorption. An iterative fit determined the continuum for each 
object. The initial input PAH template was from the NGC 7714 IRS spe ctrum after sub- 
tractin g its power-law continuum. The silicate absorption profile was from IChiar &: Tielens 
( 120061 ) with central optical depth rg 7 a free parameter. The 7.7 and 11.3 /xrn PAH line 
luminosities and equi valent widths for the local starburst sample ( IBrandl et al.ll2006l ). the 
local ULIRG s ample (lArmus et al.ll2007l ) , and the present sample were derived the same way. 
Brandl et al.l (120061 ) used a different method to derive the same parameters; their method 
would give lower 7.7/im PAH flux densities and luminosities. This is due to the complicated 
continuum at ~ 8fim. Our 11.3/im PAH flux densities are consistent with theirs. Table 3 
gives the results. 



4. Multi- Wavelength Studies of ULIRGs at 



1.9 



AEGIS (IDavis et all 120071 ) and FIDEL (IDickinson et al.l 120071 ) provide a rich X-ray to 
radio data set to study the ULIRG SEDs. Obje cts in our sample ar e measured at many key 



bands: all are detected in all four IRAC bands (iBarmbv et a. 



at 15 (llm et al.ll2008l ). and all at 1.4 GHz (llvison et al. 



20081) . all but two by AKARI 



20071 ). Most are also detected 



at 70 and 160 /im in the FIDEL survey (IDickinson et al.ll2007l ). Only two objects, EGS14 
and EGS_b2, are detected in the Chandra 200 ks X-ray imaging (iLaird et al.ll2008l ). The 
flux densities in these key bands trace stellar mass, star formation rate, and AGN activity. 



^Lii,3jj,jn and Lu^m are the 7.7 and 11.3/im PAH emission luminosities defined as Lpah = 
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Objects in the pres ent sample were also observed with MAMBO on IRAM, and most were 
detected at 1.2 mm (lYounger et al.ll2008l ). Table 4 gives the photometry, and the UV-to-radio 
SEDs are shown in Figure [71 

The multi-wavelength photometry permits us to compare the present sample with the 
sub-millimeter gal axy population. There is a small region covered by SCUBA in EGS by 
Webb et al.l (120031 ). but no galaxies in the present sample are in the SCUBA region. We fit 
FIR SEDs for the sample and predict their 850 /im flux densities Fgso to be in the range 
2.2 < Fgso < 8.4m Jy (Table 4). These values are sirnilar to the flux densities for SMGs at 
the same redshifts f lChapman et al.ll2005l : iPope et al.ll2006l . |2008| ). The median Fgso for this 
sample is 4.5 mJy, compared w ith the median Ffj^n of 5.5 raJy for SMG s at 1.5 < z < 2.2 



found by I Chapman et al.l (120051 ) and 7.5 mJy by lPope et al.l (l2006l . l2008l ). In more detail, 7 
out 12 objects in t he present samp l e hav e -Fs5 n fainter tha n 5 niJy, wh ile the flux densities 
for most SMGs in IChapman et al.l (120051 ) and iPope et al.l (l2006l . |2008| ) are brighter than 5 



mJy. We therefore argue that this sample is part of a slightly faint SMG population. 

Optical and radio morphologies of the gala xies provide impo rtant information on their 
assembly histories. HST ACS F814W imaging (ILotz et al.ll2008l ) covers the central l°xlO' 
region of the EGS. EGS l/4/b2 are outside the ACS coverage, but rough optical morpholo- 
gies are available from Subaru i?-band images. Optical images of each object are presented 
with their SEDs in Figure [71 Most objects have irregular or clumpy morphologies in the 
rest-frame 2000 A < A < 3000 A bands with a typical size of 1'.'5, suggesting extended star 
formation in a region with a size of about 13 k pc. The 1.4 GHz radio imaging of EGS has 
a mean circular beam width of ~3'.'8 FWHM (llvison et al.l 120071 ) and is unable to resolve 
morphologies except in a few cases. EGS 23 and 24 show elongated radio morphologies 
aligned with their optical extent, indicating that the radio and rest-frame UV light are from 
the same extended star formation regions in both cases. 

The spatial distributi o n of t h e stellar population i s trac ed by the rest-frame o ptical 



imagmg. [Windhorst et al.l (120021 ). iPapovich et al.l (120051 ) . and IConselice et al.l (120051 ) have 
argued that UV-dominated star-forming galaxies at high redshifts have similar morphologies 
in the rest-frame UV and optical bands. One outstanding property of galaxies in the present 
sample is their extremely red optical-NIR color. Seven objects in the sample have observed 
{R — K)Yega. > 5, qualifying them as Extremely Red Objects (ERO). EGS4 is the reddest 
with {R — K)Yega = 6.8. Red colors like these are common among distant ULIR Gs; examples 
inclu de ERO J164502+4626.4 
2002h and CFRS 14.115 7 at z 



counterparts to SMGs ( Small et al. 19991 : 



[HR9 4] 10 or sometimes "H R 10") at z = 1.44 (lElbaz et al. 
1.1 5 (iLe Floc'h et al.l 120071) . EROs are commonly seen as 



Frayer et al.ll2004l ). The red optical-NIR colors. 



corresponding to rest NUV — R for our sample, indicate either dust extinction in these 
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objects or high stellar mass. The stellar population modeling in the next paragraph suggests 
objects in our sample have both heavy dust extinction and high stellar masses. The heavy 
dust extinction does not seem to reconcile with objects being detected in the ACS 606W 
and 814W bands, which are the rest-frame 1800-2600Afor galaxies at 2; ~ 2. The irregular 
and clumpy morphologies in Figure [7| highly non-uniform dust extinction in the objects in 
our sample. Only two objects are undetected in the deep ACS 814W image, probably due 
to higher column density of dust in the compact stellar and gas distribution. 



4.1. Stellar Population and Mass in ULIRGs 



Stellar population modeling provides a way of determining physical parameters from 
the observational da ta, but it is very difficult to model stellar populations in ULIRGs. 
Tacconi et al.l (120021 ) measured dynamical masses for a sample of local ULIRGs with NIR 
spectroscopy and found stellar masses in the range of 3 x lO^^M© < < 2.4 x IO^^Mq 
with a mean of 1.1 x 10^^ Mq. Their local sample has a mean absolute K-band magnitude of 
Mk = —25.8 after adopting a dust correction of Ak = 0.7 mag. The IRAC 8 /im flux densi- 
ties of this sample correspond to a similar K-band magnitude range with < Mk >= —25.7 
if the same dust correction is used. This suggests a similar mass range for our sample, 

ULIRGs have a burst star formation history, very young stellar populations, and non- 
uniform dust distribution, all of which can introduce large uncertainties in modeling their 
stellar populations. On the other hand, stellar masses are the most robust property against 
variations i n star formation history , met allicities, and extinction law in modeling stellar 
population (IForster Schreiber et al.l |2004| ). We perform a stellar population analysis on 
the present sample, mai r ily to measure their stellar masses. We flt galaxy SEDs with 
Bruzual fc Chariot et al.l (l2003l . hereafter BC03) stella r population rnodels with a Salpeter 
IMF and a constant star forniation rate. Seve ral groups ( Shapley et al .112001 : van Dokkum et al 
2OO4J : iRigopoulou et al.ll2006l : iLai et al.ll2007l ) have argued that a constant star formation rate 
provides a reasonable description of stellar population evolution for galaxies with ongoing 
star formation at high redshifts, such as LBGs, Lyman-alpha emitters (LAEs), and DRGs. 
The stellar population age, dust reddening E{B — V), stellar mass, and derived star forma- 
tion rate from the model fltting are listed in Table 5, and the model SED flts are shown in 
Figure[71 Objects in this sample have estimated stel lar masses with > 10^^ Mp^ , similar to 



values found for local UL IRGs (ITacconi et al.ll2002l ). DRGs, and feiC galaxies (iLabbe et al. 
2005l : lDaddi et al.ll2007ah . 
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4.2. Total Infrared Luminosity and Star Formation Rate 



( Condon 


1992; 


Kennicutt 


1998) 



The validity of these methods needs to be established at high redshift. Most objects in the 
present sample are detected at 7 /xm, 160 /xm, and 1.2 mm, permitting a direct measurement 



of total infrared l uminosity Ltr (Papoyich et al.ll2007h . In practice, we derived Lm by fitting 



SED templates (jChary fc ElbazI l200ll ) to the observed 70 fim, 160 /im, and 1.2 mm flux 
densities (Fig. [7]). All galaxies in the sample have Lj^ > 10^^ Lq (Table 4), qualifying 
them as ULIRGs. EGS14 and EGS21 have Ljr > Lr^ and a re thus HyperLIRGs. 



All sample galaxies are also detected at 1.4 GHz (jivison et al.l 120071 ). We will be able to 
verify: 1) whether Liaghz is correlated with Lm for ULIRGs at z ^ 2; and 2 ) whether 



such a correlation at high redshifts is consistent with the local one JCondonlll99"i ). Figure [8] 



plots the radio luminosity Liaghz vs Lir for this sample and variety of local starburst 
and ULIRG sample s. The FIR-radio r atio q for this sampl^ are in Table 4 with a mean 



< q >= 2. 19 ±0.20. iKovacs et al.l (120061 ) measured Ljji using 350/im 850/im and 1.2mm flux 
densities and obtained a mean < q >= 2.07±0.3 for SMGs at 1 < 2; < 3. Both measurements 



yields q for ULIRGs at 2; ~ 2 close to, but smaller than the local value q=2.36. ISajina et al. 



(120081 ) showed more clearly a trend in their AGN dominated sample at 2; ~ 2: sources with 
strong PAH emission have q in 1.6 < q < 2.36; while all power-law sources have q < 1.6. 
Normally, radio excess is due to non-thermal emissio n from AGNs, but galaxy merging can 
also enhance the non-thermal synchrotron radiation (jCondonlll992l ). Merging processes are 
evident in our sample. We will argue in following paragraphs that AGNs activities may exist 
in most objects in the sample. In fact, two X-ray sources, EGS14 and EGS_b2, and the 
serendipitous power-law source EGS24a show higher radio excess (lower q) than rest objects 
in the sample. Two scenarios can be differentiated by their radio morphologies: AGNs 
are point sources and mergers, in most cases, are extended sources. Currently we cannot 
determine which scenario is responsible to the radio excess due to lower resolution of the 
1.4gHz radio images (Figure [T]). 

Another measure used to estimate Ljji for local galaxies is the IRAC 8 /im luminosity, 
Lsum, though there is considerable debate about how reliable this method is. Ls ;j,m is deflned 



as Lsu m = 47r(i|(z/F,^)8^m where F^, is the rest frame IRAC 8 /xm flux densitv (|Huang et al. 



2007bl ). Ls^m is found to be correlated with Lm for local galaxies (jWu et al.ll2005l ). The 
MIPS 24 /im band directly measures the rest IRAC 8 /xm flux densities for our sample. A 
galaxy's 8 /im flux density actually has two components (aside from starlight, which can 



log ( 3.75 l^IVm-^. 



log( 



Wm 



r) defined bv lCondonl (Il992h 
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be subtracted if necessary): the 7.7 fim PAH emission feature complex and a featureless 
continuum, coming from an AGN or warm dust in the interstellar medi um. There a re 



several models for the IR e mission from galaxies, which convert Ls/xm t o Lrn (jCharv &: Elbaz 



2001 



Bavouzet et al 



Dale fc Heloul l2002l . hereafter CEOl and DH02). Empirically, IWu et all (120051 ) and 



( 20081 ) found a correlation between Lgyum and both Liaghz and Ljr for star- 
forming galaxies. At the high luminosity end, local ULIR Gs deviate from th is correlation 
with higher i^8/im/-^^/ij ratios, such a trend was also see by iRigby et al.l (120081 ). 

Figure [9] shows a correlation between Ls^m and Lm for all populations. However, 
the Ls^m — Lir relation for objects in our sample and the local ULIRGs with high Lir 
has a higher offset than that for the loc al starburst galaxies and the model prediction 



( IChary &: Elbaa l200ll : iDale &: Heloul |2002| ) . This indicates that, for a given L/^, Lg^rn 



for objects in our sample and some of local ULIRGs is higher than the model prediction. 
Thus objects in our sample have an 8/im excess co mparing with the CEO l and DH02 model 
prediction. The empirical Ls/xm — -^/ij relation of iBavouzet et al.l (120081 ) derived with sam- 
ples at various redshifts matches local starburst galaxies, but predicts much high Ls^m for 
ULIRGs and HyperLIRGs. The Lsfxm — Ljr relation for our sample permits to estimate Lir 
for same type of objects with only 24/im flux densities. 

Our IRS spectra can be used to separate the PAH from continuum in the (rest) 8 /im 
band, and each component's contribution to Lg^m can be measured. PAH luminosity is 
thought to be a generally good tracer of star formation rate, but the Lj 7j,m / Lju ratio is 



known to be luminosity-depende nt, decreasing at high luminosity (IRigopoulou et al.lll999 



Desai et al. 


2007; 


Shi et al. 


2007) 



20071 ). Figure [TOl shows Lv vjim / Lrn versus Ltr. In this diagram. 



each population is well separated from the others. The average Lunm/ Ljr rati o for local 



ULIR Gs is seen to be lower than for local starburst galaxies. The HyperLIRGs in lYan et al. 



( 120051 ) and ISajina et al.l (l2007al ) have the lowest Lj^jj^m/ Ljr ratio. In contrast, the present 
sample has the highest Lpah/ Ljr ratio, and the trend is the same for the 11.3 /im PAH 
feature (Fig. [TTj) . Objects with such a high PAH l uminosity have neither been found locally 
nor in the MIPS 24/im luminous sample at z ~ 2 (lYan et al.ll2005l : iHouck et al.ll2005l ). 



Starburst galaxies were expected to have the highest Lrah/ Ljr, and Lpah/Lir was 
seen to decrease with increasing Ljr. Our sample shows a new ULIRG population with 
much higher PAH emissions at 7.7 and 11.3 /im. We argue that the high Lpah/ Lir ratio 
for our sample is generally compatible to extrapolation from the Lrah/ Ljr — Ljr relation 
for starburst galaxies. Both Lt j^i^ and Lii.s/im for local starburst galaxies are strongly 
correlated with Ljr in Fig. (TD] and Fig. [TTJ We flt both data sets and obtain the fol- 
lowing relations: Lir oc (i^7.7/im)°'^^ and Lir oc (i^ii.s/im)^ ''^- Both relations convert to 
Li.inm/LiR cx {LiRf-'^^ and Lii,^^xn./LiR cx {Lir)^-^^ respectively, as plotted in Fig. [Hand 
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Fig. [TTJ The Lpah-Lir relation for local starbursts predicts a higher Lpah-Lir ratio in the 
LiR range for our sample. Our sample have high Lpah-Lir ratio close to the extrapolation 
comparing with other ULIRGs population, indicating a starburst domination. The deficient 
PAH emission in our sample implies most likely existence of AGN in our sample, though 
strong UV from intensive star forming region can also destroy PAH. 

The MIR spectral properties and Lrah/Lir of our sample are closer to local s tarburst 



galaxies, even though their Ljr differs by 2 orders of magnitude. iFarrah et al.l (120071 ) reached 
the same conclusion by comparing silicate absorption strength for their sample with those for 
local ULIRG and starburst galaxies, and they propose six possible scenarios to explain the 
similarity between high redshift ULIRGs and local starburst galaxies. Our multi-wavelength 
data set provides further constrain on physical properties of our sample. The ACS I-band 
images (Figure [7]) show multi-clumpy morphologies extended to > 10 Kpc size for most 
objects in our sample. At 2; ~ 2, the observed I-band probes the rest NUV band, thus is 
sensitive to star formation. Local ULIRGs, however, have much more compact morphologies 
in the GALEX NUV images §. The extend morphologi es of our sample sup port both gas-rich 



merging and starburst geometry scenarios proposed by lFarrah et al.l (120071 ). In this scenario, 
the silicate dust column density is reduced after star formation region is stretched to a large 
scale during merging. The extended morphologies in rest NUV indicate an extended star 
formation in our sample, thus extended distribution of PAH emission. In such an extended 
distribution, more PAH can survive in strong UV radiation field from central AGN than 
those in a compact distribution. This scenario thus explains the higher Lrah/Lir in our 
sample than local ULIRGs. 

Star forming galaxies at z ^ 2 are found to generally have much less dust extinction 



than their local counterparts. iReddy et al.l (120061 ) found that there is a correlation between 



Ljr/ LiQQQ and //1600 + Ljr for star forming galaxies at 2 ~ 2, where Ligoo is the monolithic 
luminosity at I6OOA. This correlation has a higher offset than the local relation, indicating 
less dust extinction in the line of sight for galaxies at z ~ 2. Most objects ii i our sample 



l ie up on the Liqqq/ Ljr-Liqqq + Ljr relation for galaxies 2; ~ 2 (Figure [12]). iReddy et al. 



(I2OO6I ) argued that dust distribution and star formation region become more compact in local 
galaxies. We argue that lower surface density of dust density and extended star formation 
region with high SFR permit to detect both UV and PAH emission from most objects in 
ours sample. 

The star formation rate for a galaxy can be estimated from its FIR and ultraviolet 



^The GALEX UV morphologies for local starburst galaxies and local ULIRGs are from 



http: / /galex. stsci.edu/Galex View/. 
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( Kennicutt 


1998; 


Bell et al. 


2005) 



SFR/{Mq yr- 



Cx (L,^ + 3.3L28o)/i^e 



©5 



(3) 



where L2 8n is the monoch romatic luminosity (uncorrected for dust extinction) at rest frame 
280 n m JWolf et all 120051 ) . The constant C is 1.8 x 10"^° fo r the Salpeter IMF jKennicutt 



19981 1 ■ and 9.8 x 10"^^ for the Kroupa IMF flBell et all 120051 1 . In the following text, we will 
adopt the Salpeter IMF for all L^/j-SFR conve rsion in this pa. per. SFR will reduce by a 
factor of ~2 if we switch to the Kroupa IMF (IBell et al.l 120051 ) . The 280 nm band shifts 
to the observed / band at z ~ 2. L280 was calculated from the ACS F814W magnitude if 
available or otherwise the CFHT / magnitude. All objects in our sample have L280 in the 
range 5 x 10^ < L280 < 3 x lO^^L©, less than 1% of their Lj^. The star formation rate seen 
at rest-frame 280 nm is at most 20 Mq yr~^, and most UV light produced by newborn stars 
is absorbed by dust and re-emitted in the infrared. Thus we omit the L280 contribution in 
our SFR calculation. 



Total infrared luminosity, of ULIRGs may be partly powered by AGNs (INardini et al 



20081 ) . thus using L/^j may over-estimate their SFR. The PAH emission only traces star for- 
mation, and is free of AGN contamination. We calculate SFR for our sample with their 
LpAH using the Lpah — SFR relation, established from local starburst galaxies shown in 
Figure [To] and Figure [TTl Results are given in Table 5. Star formation rates for our sample 
converted from Lm using Equation 3 are much higher, with an average SFR ~ 1000 Mq yr^^. 
-^7.7/im and Ivii.3^rn (Table 5) give smaller, star formation rates, in the range 150 < SFR < 
600 Mq yr~^ for most objects, that are quite consistent with the stellar population modeling 
results. The discrepancy between both star formation estimations may be due to: 1. part of 
star formation occurs in region with no PAH, thus Lpah underestimates the SFR; 2. Ljp 
contains AGN contribution, thus over-estimate the SFR. It is very possible that both can 
happen in one object simultaneously, namely its AGN destroys PAH in surrounding area 
where star formation occurs. This will further increase the discrepancy, so the real SFR 
should be in between both estimations. 

Our sample have both high star formation rate and high stellar mass, supporting the 
galaxy formation in the "downsizing" mode. The star formation rates and stellar masses for 
our sample are cor isistent with the SF R-stellar mass relation obtained from BzK galaxies at 
z ^ 2 (Figure [T3!l. iDaddi et al.l (l2007al ) sh owed that simulated ga l axy po pulat ions taken from 



the M ilUennium Simulation lightcones of iKitzbichler fc White I (120071 ) and ICattaneo et al. 
(120081 ) failed to re-produce the SFR-star mass relation at z = 2, thus underestimate number 
of ULIRGs at 2 ~ 2. 



It has been long anticipated that ULIRGs have a dominant co ntribution to the tota l 
infrared luminosity density, thus star formation rate density, at z ~ 2 (ILe Floc'h et al.ll2005l ). 
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We use the Vmax method to calculate the total infr a red lu minosity density for our sample to 
be 2.6 X 10"^ Lq/Mpc^. The sample of iFarrah et al.l (120071 ) with the same limiting flux yields 
a density of 8.8 x IO^Lq/Mpc^. We argue that the difference is due to the cosmic variance, 
because these objects are massive galaxies and thus have a much stronger spatial correlation. 
Both densities are low er than ULIRG Lm density at 2; ~ 1, ~ 10^ Lp^/Mpc^ for all UL IRGs 
(ILe Floc'h et al.ll2005l ). Most objects in our sample and those of iFarrah et al.l (120071 ) have 
LiR > 5 X lO^^L©, the major c ontribution to the LrR density at ^ ~ 2 comes from ULIRGs 
with 10^2 <LiR<5x 10^2^0 (iHuang et al.lboooh . 



4.3. AGN in the z ~ 1.9 ULIRG sample 



One direct way of identifying an object as an AGN is to measure its X-ray luminosity. 
Two objects in the sam ple, EGS14 and EGS_b2, are in the rn ain AEGIS-X catalog from the 
Chandra 200 ks images (jNandra et al.ll2007l : iLaird et al.ll2008l ). Their X-ray fluxes Fo.5_io kev 
are 1.2 x 10~^^ and 6 .4 x 10~^^ erg cm~^ s~^ r espectively. Calculated X-ray luminosities Lx 



(INandra et al.l 120071 : iGeorgakakis et al.l l2007h are 1.0 x 10^=^ erg s-^ for EGS14 and 9.4 x 



Xq43 gj,g g-1 i^Qj, EGS_b2. Hardness ratios are 0.45 and -0.30, respectively. Therefore EGS14 is 
a type 2 (obscured) AGN, and EGS_b2 is ver y close to a type 1 ( unobscured) QSO according 
to the X-ray luminosity and hardness ratios (jSzokoly et al.ll2004j ). In addition to EGS14 and 
EGS_b2, EGSl has a low-significance X-ray counterpart. At the location of this source there 
were 6.5 net soft band counts (10 counts total with an estimated 3.5 count background). This 
gives a Poisson probability of a false detection of 3.5 x 10~^. The source was not detected 
in the hard band. If the detection is real, EGSl has -fo.5-2 keV 

with L2-IO keV 



3.0 X 10 erg cm ^ s ^ 



1.1 X 10^^ erg s ^, thus is qualified to be an AGN. 



The remaining 10 ULIRGs are not detected in the current Chandra observation. Stack- 
ing in the soft band gives 19.5 counts above an average background of 9.85, corresponding 
to -Fo.5-2 keV = 3.8 X 10^^^ erg cm~^ s~^ or Lx = 1.1 x 10^^ erg s~^ at 2a significance. There 
was no detection in the hard band. Even if EGSl is added to the stacking, nothing shows 
up in the hard band, but the soft band detection significance rises to 3.2cr. The mean fiux 



is -Fo.5-2 keV = 4.7 X lO^"*^"^ erg cm~^ s^^ or Lx = 1.3 x 10^^ erg s~^. This average X-ray 
lu minosity represents either a very weak AGN or strong star formation. Using the relation 
of iRanalli et al.l (120031 ). this average X-ray luminosity corresponds to a star formation rate 
of 220 Mq/yt, consistent with the SED and PAH estimation. However, we argue that the 
stacked X-ray signal comes from central point sources. These objects have very extended 
and elongated morphologies in the rest-frame NUV band. If the X-ray photons are from 
these star formation regions, stacking would not yield any signal unless they are aligned. 
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Emission in the rest 3-6 um wavelength range is anothe r indicator of AGN activity 



JCarleton et al.lll987 



2006 



Shi et al.ll2007l : lAlonso-Herrero et al.ll2006l : iHines et all2006l : Ijiang et al 



Shi et al.l 120071 ). The longer end of that range, which has minimal stellar and PAH 



emission contamination, is ideal for detecting what is nowadays thought to be hot dust emis- 
sion closely related to the AGN accretion disk. Luminosity at these wavelength s (L5) can 
be converted to Lj r for QSOs wi th the QSO SED templates JeIvIs et al.lll994h . AKARI 
15 ixm. photometry (jim et al.ll2008l ) provides the best measurement of L5 for our sample. All 
galaxies within the AKARI c overage are det ected except EGS26, for which the 3ct limiting 
flux density is F15 < 58 /iJy Jim et al.l[2008l l The AKARI 15 yum band is wide enough to 
include the 6.2 /im PAH feature for objects with 1.5 < z < 2.2, but this feature is much 
weaker than the 7.7 yum feature. Thus the AKARI 15 yum band is a better measure of AGN 
emission than the MIPS 24 yum band. 

In fact, the -F15/F24 ratio for our sample measures the continuum-to-PAH ratio, and 
thus the AGN fraction. Figure \\M shows this ratio versus redshift. The ratios for the two 
known AGNs with AKARI coverage, EGS14 and EGS24a, are very close to the expected 
values for Seyfert 2's. EGSll and EGS12 are similar to expectations for H ll-type ULIRGs. 
The flux ratios for the remaining objects in our sample sh ow even more PA H than starbursts, 
indicating starburst domination in these objects. SMGs (jPope et al.ll2008l ) have very similar 
F15/F24, ratios as objects in the present sample, implying the same properties shared by both 
samples. The SMGs also show very strong PAH features in their IRS spectra. This supports 
the argument that most objects in the present sample are part of a SMG population, and 
are starburst dominated ULIRGs. 

A starburst dominated ULIRG can still have a deeply dust-obscured AGN. Many cur- 



rent theoretical models (e .g. , iMihos fc Hernquistll 19941 . 1 19961 : iDubinski et al.lll999l : ICox. Jonsson. Primack. I 
20061 : iHopkins et al.ll2006l ) suggest that such a dust-obscured AGN can have a significant con- 
tribution to Ljr of a ULIRG. A study of l ocal ULIRG IRS sp e ctra show an average of 15% 
LiR are from central dust-obscured AGNs (INardini et al.ll2008l ) . iNardini et al.l (120081 ) argued 
that ULIRG luminosity in 5 < A < 6 is dominated by hot dust emission from AGNs. Most 
objects in the present sample are detected at 15/im thus permit to measure their rest 5/im 
luminosities, Ls^m? which trace AGN activity. Ls^ni for the present sample is in range of 
9.9 < Loq(Lr,/,,m/Lcr,) < 12.6 (Table 4). Using (Ljr/ L^^m)QSO = 22.8 from the IeIvIs et aL 



(I1994J ) QSO SED, we calculate that such a QSO contribution is about 14% of Ljr for objects 
in our sample, consistent with those for local ULIRGs. 
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5. Summary 

The results for the present sample combined with others in Table 1 show that high- 
redshift ULIRGs have a diverse range of properties, and different selection criteria to pick 
out different populations. The combination of IRAC colors and MIPS 24 fim flux used here 
selects ULIRGs with strong 7.7 /xm PAH in a rather narrow redshift range around z ~ 1.9. 
This sample shows a starburst dominated stage in gas-rich merging powered ULIRGs at 
2; ~ 2. In this stage, intensive star formation occurs in much extended region with a typical 
scale of ~ 15 Kpc indicated by their ACS morphologies. Objects in this sample have higher 
total infrared luminosities than local ULIRGs, but the Lpah/Ljr ratios for the sample are 
higher than those of local ULIRGs. We argue that the high Lpah/Ljr ratio is due to the 
extended PAH distribution, which is less affected by strong UV emission from central AGNs. 
Most objects follows the same Ljr/ LiQOQ-Lboi relation as that for BM/BX, DRG and BzK 
galaxies, though they are at higher luminosity end. 

Stellar masses in this sample already exceed 10^^ Mq. Most stars must have formed 
prior to this stage. The SFR-stellar-mass relation for this sample is also consistent with 
that for the rest populations at 2; ~ 2, which is much higher than the theoretical model 
prediction. 

Only a few of the ULIRGs in our sample show direct evidence to have AGNs with either 
high X-ray luminosities or hot dust emission in the mid-infrared. Several pieces of evidence 
show weak AGNs existing in this starburst dominated ULIRG sample: systematically higher 
LiAGHz/ LjR ratio than the local radio-FIR relation, and an average X-ray emission of Lx = 
1.3 X lO'*^ erg s~^ from point sources. AGN contributes on average 15% of total infrared 
luminosity for our sample. 

This sample presents an early stage with very intensive star formation but weak or heav- 
ily obscur ed AGNs. ULIRGs in other samples at similar redshift but with different selection 



methods (lYan et al.l l2005l : ISajina et al.l l2007al ) have higher total infrared luminosities and 



lower PAH luminosities, indicating increasing AGN and decreasing star formation at higher 

LiR. 



This work is based in part on observations made with the Spitzer Space Telescope, 
which is operated by the Jet Propulsion Laboratory, California Institute of Technology under 
a contract with NASA. Support for this work was provided by NASA through an award 
issued by JPL/Caltech. 
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Fig. 1. — IRAC color-color diagram for several samples. The upper left panel is for the 
entire AEGIS spectroscopic redshift sample with < 2; < 1.5; the upper right p anel is for the 



comb i ned BM/BX and L BG samples with confirmed spectroscopi c redshifts ([Steidel et al. 



2OO4J : iReddv et all 120081 ): the lower left p anel is for the DR Gs flFranx et all l2003f ): and 
the lower right panel is for BzK galaxies (jPaddi et al.l 120051 ). The boxes in each panel 
denote the IRAC color selection for the present sample. The track for the M82 template 
is also plotted in each panel. The rest-frame UV-selected BM/BX and LBG galaxies have 



generally faint IRAC flux den sities and thus larger photometric uncertainties (IHuang et al. 



2OO5I : iRigopoulou et al.ll2006l ). increasing the apparent scatter in the upper right panel. 
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Fig. 2. — IRAC color-color diagram for EGS galaxies with F{24 /im) > 80 /iJy. Small 
dots show all such galaxies; large red dots show galaxies in the current IRS spectroscopic 
sample, which also requires F{24 fim) > 500 /iJy. The black box shows the IRAC color 
criteria, which should select objects at 2; > 1.5. The one red dot outside the selection box is 
the serendipitous source EGS24a. Objects from other IRS spectroscopic samples (Table 1) 
at ^ ~ 2 are plott ed for comparison: blue triangles a nd green diamonds denote "optically 



starbursts of Weedman et al 



invisible" so urces flHouck et al.l l2005l : lYan et al.l l2005l ) . Blue squares denote the luminous 



J2006bh . 
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Fig. 3. — Observed IRS spectra. The vertical scale is linear but different for each panel. The 
gray-scale images are the two-dimensional IRS spectral images after wavelength and position 
calibration. Each image shows 5 pixels or 25'.'5 along the slit. The dashed lines indicate the 
central wavelengths of the PAH emission features, rest-frame 7.7, 8.6, and 11.3 /xm left to 
right. EGS24a, the serendipitous object in the slit of EGS24, shows a power-law SED with 
strong silicate absorption. Cross-correlation of the template (red dashed line) of the local 
ULIRG IRAS F08572+3915 to the spectrum of EGS24a gives z = 2.12. For EGS_b2 at 
z = 1.59, the 7.7 fim feature is off scale to the left, and the peak observed at 33.5 fim is the 
[Ne II] emission line at rest wavelength 12.81 /xm. EGS_b6 is a confused case with combining 
spectra of two galaxies at z=1.02 and z=2.0. 
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Fig. 4. — Near infrared spectrum of EGS12 taken with the MOIRC spectragraph on Subaru. 
There is nothing detected except one emission hne at 1.992 /im. We identiiy this hne as Ha 
dX z = 2.033, and corresponding rest-frame wavelengths are marked above the plot. The 
redshift is consistent with z = 2.03 derived from the PAH features in the IRS spectrum 

(Fig. ED. 
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Fig. 5. — Redshift distributions for spectroscopic samples in Table 1. In the third panel, the 
red line shows the distribution for the starburst (SB) sample and the black line for the AGN 
sample. 
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Fig. 6. — Stacked spectrum for ULIRGs in the present sample (black line). The short 
wavelength limit for the stacked spectrum is 7 /xm. The large dot at rest wavelength 5.3 /xm 
represents the stacked AKARI 15 /xm flux density. The vertical scale is linear flux density per 
unit frequency in arbitrary units. Other lines show s tacked spectra of comparison samples: 
local starburst galaxies (grey line: lBrandl et al.ll2006l ). local Seyfert-type ULIRGs (red line), 
local LINER-type ULIRGs (blue line), and local H ll/sta rburst ULIRGs (green line). The 



local ULIRG samples are from the IRAS 1 Jy sa mple flKim &: Sanders.l 



observations in the IRS GTO program (PID 105: IParrah et al.l booi. krmus et al.l boOTh. 



1998h with IR S 



Types were assigned according to optical spectroscopy (jVeilleux et aL 1995 . 1999 ). The 



average SED of the present sample is ver y similar to those of local LINER and H ll-type 
ULIRGs, while average SEDs for objects in lSajina et al.l (j2007al ) are close to the local Seyfert 
tvDe ULIRGs with much hieher continuum emission. 
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Fig. 7. — Spectral energy distributions and morphologies for sample galaxies. Inset images 
are in negative grey scale and are 12" square. The red images come from HST ACS data 
(filter F814W) if available; otherwise Subaru R. Black dots represent photometric data, and 
the blue line is the stellar population model (BC03) that best fits each source. The red lines 
are the CEOl dust templates chosen to match the FIR luminosity of each source. 
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Fig. 8. — Correlation between Ljr and Li^ghz- Ljr was calculated by fitting SED templates 
(CEOl) to the MIPS 70, 160 /zm, and MAMBO 1.2 mm photometry. The two points with 
red circles are the X-ray sources EGS14 and EGS_b2, and the one with a blue circle is the 
serendipitous object EGS24a. The sample is plotted together with local starburst galaxies 
and ULIRGs against the the local FIR- Radio relation (the thick line, Condon, 1992). The 
color coding: Seyferts in the local ULIRG and starburst samples are shown in red; LINERs 
in the local ULIRG and starburst samples are shown in green; and Starburst/HII-type 
ULIRGs in blue, the same as in Figure [61 This plot shows that the local starburst galaxies 
and ULIRGs, and objects i n our sample are all consistent with the local FIR- Radio relation 
(Condon, 1992). Objects in lSajina et al.l (j2007al ) show strong radio excesses indicating AGNs 
in their sample. 
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Fig. 9. — Rest-frame 8 /iin luminosity, Lg /im versus total infrared luminosity, Ljji. ULIRGs 
in this sample are shown as filled black dots; the two points with red circles are the X- 
ray sources EGS14 and EGS_b2, and the one with a blue circle is the serendipitous object 
EGS24a: local starburst galaxies as triangles, and local ULIRGs as diamonds. Local star- 
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Fig. 10. — The 7.7 /iin PAH to total infrared luminosity ratio versus total infrared luminosity. 
The L7 7/L7/J ratio measure the star formation contribution in the total infrared luminosity 
for objects in the sample. The present sample has the highest L-j -j/ Ljji ratio, indicating that 
they are starburst dominated ULIRGs. The L-j^t/Lhi ratio for the present sample is still 
compatible to the empirical relation of Ly.y/L//? ~ LPf^ from the local starburst galaxies. 
Objects in the present sample sample are shown as filled black dots (the X-ray source EGS14 
is a filled dots with red circle, the 7.7/im PAH in the other X-ray sources EGS_b2 is not in 
our observation band, thus it is not in the diagram); local starburst galaxies as triangles, and 
local ULIRGs as dia monds. Local starbu r sts and ULIRGs are colo r-coded based on their 
spectral classification f Veilleux et al. IQQsl . 1999 : Brandl et al. 2006 ): red for Seyfert, green 
for LINER, and blue for H ll/starburst. Ope n circles show ULIR Gs and HyperLIRGs from 
the lYan et al.l ( l2005l ) sample with data from ISajina et al.l (l2007al ). The inserted plot shows 
a strong correlation between L7 7 and for the local starburst galaxies. Thick lines show 
linear fits to the correlations for local starburst galaxies as Lm ~ -^^77^, which transfers to 
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Fig. 11. — The 11.3 /im PAH to total infrared luminosity ratio versus total infrared lumi- 
nosity. The plot shows the same pattern as in Figure [TOl Objects in the present sample 
sample are shown as filled black dots; the two filled dots with red circles are two X-ray 
sources, EGS14 and EGS_b2; local starburst galaxies as triangles, and local ULIRGs as di- 
amonds. Local starbursts and ULIRGs are c olor-coded based their spectral classification 
(IVeilleux et aDllQQsl . Il999l : ISrandl et allboOGh : red for Seyfert, green for LINER, and blue 
for H ll/starburst. The inserted plot shows a strong correlation between L11.3 and Lm for 
the local starburst galaxies. Thick lines show linear fits to the correlations for local starburst 
galaxies as L/^j ~ Llf^^, which transfers to Lh^^/Lir ~ L^j'^ plotted with the thick line. 
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Fig. 12. — The IR-to-UV luminosity ration for galaxies at 2; ~ 2 and z ^ 0. The left panel 
is for galaxies at z ~ 2. Our sample is plotted against galaxies at z ~ 2 selected in various 
bands. The solid lines in both panels are the Ljul L^mn-Lju + Ligoo relation for BM/BX 
sources, fei^ galaxies, DRGs, and SMGs (IReddv et al.ll2006l ). The rig ht panel is for loca l 
galaxies including normal galaxies (IBell et al.ll2005l ). starburst galaxies (IBrandl et al.ll2006l ). 
and ULIRGs. Most objects in our sample have the same relation as the rest of galaxies 
population at 2: ~ 2. Three obje cts in our sample with extreme red colors, together with 
some DRGs in lReddy et al.l (120061 ) . are off the relation. They locate in the region where local 
ULIRGs are, indicating a compact dust distribution in those objects. See detailed discussion 
in the text. 
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Fig. 13. — The SFR-stellar mass relation for galaxies at z=0, 1, and 2 suggests the "down- 
sizing" scenario for galaxy formation. The mean Lm-M^ relation for BzK at 2; ~ 2 is from 
Daddi et al.l (12007ar). the relations for GOODS galaxies at 2; ~ 1 and SDSS galaxies at 



2; ~ 0.1 are from lElbaz et al.l (120071 ). Again, objects in our sample are consistent with the 
relation for B M/BX sources and BzK , but a t hig h mass end. Bo t h simu lated galaxy popula- 
tion models of iKitzbichler fc White I (120071 ) and ICattaneo et al.l ( l2008l ) predicts much lower 
star formation rate for galaxies with a given stellar mass. 
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Fig. 14. — F(15 Mm) 7^(24 M m') versus redshifts for the ULIRG sample (filled circles). The 
SMGs from jPope et al.l (120081 ) are also plotted (red open circles). This ratio should measure 
the continuum (hot dust) to PAH ratio for objects in this redshift ranger and hence the AGN 
contribution to the mid-IR luminosity. Lines show the relations defined by local templates 
(Fig. [6]): starburst galaxies (grey), Seyfert-type ULIRGs (red), LINER-type ULIRGs (green), 
H ll/starburst ULIRGs (blue), PG QSO (dashed line). One X-ray source, EGS_b2, is outside 
the AKARI15 fim coverage and not plotted. The other X-ray source, EGS14 (red diamond), 
and the serendipitous object, EGS24a (blue triangle), have colors consistent with AGN. The 
16 yum flux densities for SMGs are measured from the IRS peak-up imaging. Since the IRS 
16 yum peak-up filter profile is very similar to the AKARI 15 /xm filter profile, no correction 
is applied to the IRS 16 fim flux densities for SMGs. 
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Table 1 
IRS Observation Parameters 



Nickname^ 


EGSIRAC^ 






RA 


Dec 

J2000 




F(24/im) 

niJv 
111 J J 


Cycles 


Exp time 
s 


EGSl 


J142301. 49+533222 


4 


14 


23 


m 


50 


+53-32 


22 


6 


0.55 


10 


7314 


EGS4 


J142148.49+531534 


5 


14 


21 


48 


49 


+53:15 


34 


5 


0.56 


10 


7314 


EGSIO 


J141928.10+524342 


1 


14 


19 


28 


09 


+52:43 


42 


2 


0.62 


8 


5851 


EGSll 


J141920.44+525037 


7 


14 


19 


17 


44 


+52:49 


21 


5 


0.59 


8 


5851 


EGS12 


J141917.45+524921 


5 


14 


19 


20 


45 


+52:50 


37 


9 


0.74 


5 


3657 


EGS14 


J141900.24+524948 


3 


14 


19 


00 


27 


+52:49 


48 


1 


1.05 


3 


2194 


EGS23 


J141822.47+523937 


7 


14 


18 


22 


48 


+52:39 


37 


9 


0.67 


7 


5120 


EGS24 


J141834.58+524505 


9 


14 


18 


34 


55 


+52:45 


06 


3 


0.66 


7 


5120 


EGS24a= 


J141836.77+524603 


9 


14 


18 


36 


77 


+52:46 


03 


9 


0.66 


7 


5120 


EGS26 


J141746.22+523322 


2 


14 


17 


46 


22 


+52:33 


22 


4 


0.49 


11 


8045 


EGS_b2 


J142219.81+531950 


3 


14 


22 


19 


80 


+53:19 


50 


4 


0.62 


8 


5851 


EGS_b6 


J142102.68+530224 


5 


14 


21 


02 


67 


+53:02 


24 


8 


0.72 


6 


4388 



"Nicknames are the target names in the Spitzer archive and are used for convenience in this paper, but they 
are not official names and should not be used as standalone source identifications. 
''Source name from Barmbv et al. ( 2008h . 

^Serendipitous source found in the slit while observing EGS24. 



Note. — RA and Dec are the commanded telescope pointing coordinates. Telescope pointing was based 
on high accuracy peakup on nearby 2MASS catalog objects with the blue peakup array. Ramp duration was 
120 s for all objects. 



Table 2 

IRS Sample Selection Criteria 



Sample 24 jjm flux density Color criteria 



Houck et al. f2005) 


>0.75 


mJy 


z/F^(24/im)/i/F^(/) > 60 


Yan et al. r2nn5) 


>0.9 


mJy 


uF^{2Afim)/uF^{I) > 10 and 








uF^{24fim)/uF^{8fim) > 3.16 


Weedman et al. r2nn6bVAGN) 


>1.0 


mJy 


F(X-ray)'' > 10~^^ erg cm"^ s"^ 


Weedman et al. f2006bVSB) 


>1.0 


mJy 


IRAC flux density peak at either 4.5 or 5.8/im 


This paper 


>0.5 


mJy 


< [3.6] - [4.5] < 0.4 and 








-0.7 < [3.6] - [8.0] < 0.5 



Chandra 0.3-8 keV flux density 
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Table 3 

pah properties for the sample 



Object 


redshift" 


redshift'' 


loeL('7 7) 


EWiy.7) 


logLfll 3) 


EW{11.3) 
















Ln 








/xm 




EGSl 


1 


.95±0, 


.03 


1 


.90±0, 


.02 


11 23±0 03 


2 38±0 22 


10 18±0 15 


1. 


,68±0, 


.26 


EGS4 


1, 


.94±0, 


.03 


1 


.88±0, 


.02 


10 89±0 06 


57t0 07 


9 82±0 37 


0. 


,17±0, 


.07 




1, 


.94±0, 


.02 


1 


.94±0, 


.01 


1 1 QQ_i_n no 


o QO-un 1 
z.ov±\j.l.z 


1 n n/i -Un qi 


0. 


,26±0, 


.10 


EGSll 


1 


.80±0, 


.02 


1 


.80±0, 


.01 


11.02±0.05 


0.79±0.10 


10.25±0.12 


1. 


,19±0, 


.16 


EGS12 


2, 


.01±0 


,03 


2, 


.02±0, 


,03 


11.37±0.02 


1.46±0.08 


10.61±0.11 


1. 


,28±0, 


.55 


EGS14 


1, 


.87±0, 


.06 


1, 


.86±0, 


.03 


11.33±0.04 


1.13±0.09 


10.63±0.10 


2. 


,98±0, 


.35 


EGS21 


3 


.01±0, 


.03 


3 


.00±0, 


,03 


11.73±0.06 


1.59±0.10 










EGS23 


1, 


.77±0, 


.02 


1, 


.77±0, 


,01 


11.15±0.04 


1.45±0.12 


10.54±0.05 


1. 


,08±0, 


.08 


EGS24 


1 


.85±0, 


.03 


1 


.85±0, 


.01 


11.25±0.03 


2.24±0.18 


10.56±0.07 


0. 


,36±0, 


.08 


EGS26 


1, 


.77±0, 


.03 


1 


.78±0, 


.02 


11.16±0.03 


2.61±0.20 


10.42±0.06 


1. 


,12±0, 


.18 


EGS_b2 


1 


.59±0, 


.01 


1 


.60±0, 


.01 






10.45±0.04 


0. 


,30±0, 


.04 



Note. — • • • indicates the PAH feature lies outside the observed spectral coverage. EGS24a shows 
silicate absorption and no PAH omission (Fig. 3). 

"redshifts obtained with a ULIRG template. 

''redshifts obtained with a starburst template. 
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Tablb 4 

IR/Radio flux and luminosity of the Sample 



Object F(3.6 Atm) F(4.5 Atm) F(5.8 ^m) F(8.0 fim) F(15 urn) F(24 ^tm) F{70 fim) F(160 ^jm) F(850 fim)'=F{l.l mm) F(1A GHz) 









MJy 


MJy 


MJy 


MJy 


m.Ty 


mJy 


mJy 


mJy 


mJy Lq 


EGSl 


45.0±0.3 


55.4±0.4 


63.6±1.5 


56.3±1.6 


. . . d 


554±35 


<1.5 


12.1±8.9 


3.3 


1.86±0.50 


0.069±0.010 12.72±0.15 


EGS4 


32.1±0.3 


44.9±0.4 


52.1±1.5 


40.1±1.5 


125±24 


557±22 


2.4±0.5 


<21.0 


3.9 


1.87±0.48 


0.062±0.010 12.62±0.12 


EGSIO 


21.8±0.3 


23.0±0.3 


34.0±1.4 


28.9±1.5 


77±28 


623±35 


4.2±0.7 


45.5±8.7 


5.2 


1.65±0.69 


0.085±0.014 12.83±0.10 


EGSll 


27.8±0.3 


36.9±0.4 


38.2±1.4 


30.6±1.5 


192±31 


591±20 


5.0±0.6 


<21.0 


3.3 


0.85±0.44 


0.067±0.017 12.58±0.09 


EGS12 


27.7±0.3 


31.0±0.3 


38.3±1.4 


33.2±1.5 


196±25 


743±23 


3.9±0.6 


a 


5.4 


1.58±0.47 


0.036±0.010 12.77±0.07 


EGS14 


66.1±0.2 


89.6±0.4 


101.7±1.5 


88.4±1.6 


457±39 


1053±41 


3.8±0.6 


76.7±9.6 


6.4 


4.54±0.68 


0.316±0.023 13.18±0.06 


EGS21 


39.5±0.3 


45.3±0.4 


50.7±1.5 


35. Oil. 5 


.59±14 


605±23 


2.8±0.5 


34.2±9.4 


8.4 


1.31±0.35 


0.070±0.014 13.15±0.07 


EGS23 


47.8±0.3 


60.6±0.4 


69.1±1.5 


51.3±1.5 


132±29 


665±18 


3.7±0.4 


62.4±8.7 


4.5 


1.81±0.40 


0.119±0.015 12.79±0.08 


EGS24 


36.4±0.3 


44.0±0.4 


46.8±1.4 


37.1±1.5 


65±25 


663±29 


3.4±0.6 


9.7±9.0 


2.7 


1.49±0.74 


0.047±0.012 12.51±0.18 


EGS26 


31.7±0.3 


43.3±0.4 


47.1±1.4 


34.0±1.5 


58±20 


492±16 


1.5±0.5 


21.6±8.4 


4.5 


1.14±0.36 


0.097±0.017 12.49±0.15 


EGS24a 


22.3±0.3 


32.3±0.3 


46.7±1.5 


575±1.6 


223±36 


997±30 


2.5±0.5 


15.1±8.4 


6.4 


2.87±0.54 


0.112±0.013 12.91±0.10 


EGS.b2 


94.0±0.2 


124.8±0.4 


115.0±1.5 


117.1±1.6 




616±30 


3.4±0.5 


21.7±7.0 


2.2 




0.151±0.009 12.34±0.14 



"Confused 

''Prediction from the SED fitting. 

^FIR luminosity of the best-fit CEOS template. 

''not observed 
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Table 5 

Stellar Population Fitting Parameters 



Name 


Aere 


E{B - V) 




SFR 




Gvr 




10^^ Mq 




EGSl 


1.9 


0.3 


5 


240 


EGS4 


1.4 


0.7 


5 


320 


EGSIO 


1.1 


0.4 


2 


182 


EGSll 


2.0 


0.7 


4 


196 


EGS12 


0.29 


0.4 


1 


480 


EGS14^ 


0.26 


0.6 


3 


1320 


EGS23 


1.1 


0.6 


5 


400 


EGS24 


0.29 


0.5 


5 


580 


EGS26 


1.8 


0.6 


4 


220 


EGS_b2^ 


0.03 


0.6 


0.9 


3800 



"EGS14 and EGS_b2 are X-ray sources, and their SEDs may be contaminated by AGNs. 

Note. — EGS24a. the serendipitous object next to EGS24, was not fit with any model because it has a 
power-law SED in the IRAC bands. 
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Table 6 

Comparison of Star Formation Rates 
Name SFR(BC03)^ SFR(7.7)^ SFR(11.3)^ SFR(L/r)^ 



EGSl 


240 


386±18 


263± 68 


945±326 


EGS4 


320 


226±21 


142± 90 


750±207 


EGSIO 


182 


452±14 


207±110 


1217±280 


EGSll 


196 


277±21 


297± 61 


684±142 


EGS12 


480 


481±15 


549±103 


1060±171 


EGS14 


1320 


451±28 


568± 97 


2724±376 


EGS23 


400 


340±21 


487± 42 


1110±204 


EGS24 


580 


398±19 


504± 60 


582±241 


EGS26 


220 


346±16 


390± 41 


556±192 


EGS_b2 


3800 




417± 29 


394±127 



"Star formation rate calculated from stellar population model (BC03) fitting 
''Star Formation Rate calculated from PAH feature luminosity. 
'^Star Formation Rate calculated from far infrared luminosity Lm 




10 11 12 13 
Log(U/LQ) 



